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LOW COMPLEXITY CHANNEL ESTIMATION FOR ORTHOGONAL FREQUENCY 
DIVISION MODULATION SYSTEMS 

CROSS - REFERENCE TO RELATED APPLICATION 

[0001] This application claims priority to U.S. 
Provisional Application Serial No. 60/446,795, filed on 
February 12, 2003 and entitled, "Apparatus And Method For 
Low Complexity Channel Estimation For OFDM Systems" . 

BACKGROUND 

[0002] Fig. 1A illustrates a wireless communication 
system 13 0, which includes a first transceiver 100 and a 
second transceiver 102. The first transceiver 100 may be 
designated a "transmitter" because it first transmits 
signals to the second transceiver 102, which may be 
designated a "receiver." Both transmitter 100 and receiver 
102 may transmit and receive wireless signals, as shown by 
the transmit portions 101A, 101B and receive portions 103A, 
103B. 

[0003] Orthogonal frequency division modulation (OFDM) 
is a modulation technique for communications which splits a 
data stream into multiple radio frequency channels, which 
are each sent over a " subcarrier" frequency. 

[0004] In an OFDM baseband system, such as an IEEE 
802.11a/g system, a transmitter 100 transmits two identical 
training sequences known as "preambles" to a receiver 102. 
The transmitter 100 sends these training sequences over two 
OFDM symbol durations. The transmitter 100 may send 
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training subcarriers (also called "carriers" or "tones") 
for the training sequences that correspond to all data 
subcarriers to the receiver 102. 

[0005] Fig. IB illustrates an example of a data 
burst/packet 120 transmitted by the transmitter 100 to the 
receiver 102 in Fig. 1A. In IEEE 802.11, the packet 120 
contains a preamble 121A,B in a training mode and a 
plurality of data symbols 122A-122C in a data mode. In 
IEEE 802.16a, the packet contains a preamble 124 in a 
training mode and a plurality of data symbols 12 6A-C in a 
data mode. 

[0006] Fig. 1C illustrates a set of subcarriers for a 
preamble "P" 204 at the transmitter 100 in Fig. 1A. Fig. 
1C shows how the subcarriers are transformed and 
transmitted to the receiver 102. A subcarrier "i" in the 
frequency domain contains a training symbol . The 
transmitter 100 performs an inverse fast Fourier transform 

(IFFT) to transform the subcarriers to one time -domain OFDM 
symbol. One OFDM symbol contains multiple training symbols 
carried over the sub-carriers. The transmitter 100 
transmits the OFDM symbol across a channel. The receiver 
102 receives the OFDM symbol and performs an FFT operation. 
The receiver 102 estimates a frequency response Hi of the 
subcarrier i . 

[0007] The receiver 102 typically uses two identical 
received training sequences/preambles to compute their 
correlation and obtain time and frequency information for 
frequency synchronization and channel estimation. The 
receiver 102 may compute channel estimation using a Least 
Square (LS) estimator. An estimated channel frequency 
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response H k corresponding to a subcarrier k (k is an index 
for subcarriers) may be expressed as: 

(1) H k = ^ for k=0, 1„, N-l, and P k * 0, 

where N denotes the total number of subcarriers of an 
OFDM symbol, P* is the received symbol of the kth 
subcarrier, and P k is the corresponding training symbol for 
that kth subcarrier, which forms a part of a preamble. 
Equation (1) indicates that a channel estimate may be 
determined for each subcarrier k. 

[0008] A preamble represents information about packets 
that follow the preamble. As such, preambles are overheads 
of a packet to be transmitted. There may be multiple 
preambles. Preambles reduce the throughput of a 
communication system. For a diversity system, such as a 
multiple-in-multiple-out (MIMO) system with multiple 
transmit and receive antennas 104, 106, channels between 
all pairs of the transmit and receive antennas 104, 106 
should be estimated. Extra training overheads may be 
required. 

[0009] It may be advantageous to minimize the number of 
preambles required to be transmitted while still obtaining 
satisfactory system performance. To increase system 
throughput and reduce the overhead of transmitting two 
identical training sequences over two OFDM symbol durations 
per packet, the IEEE 802.16a OFDM system proposes 
transmitting two identical training sequences using one 
preamble over one OFDM symbol duration. The receiver 102 
uses the two identical sequences for correlation to obtain 
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frequency information. To obtain two identical training 
sequences in the time domain over just one symbol duration, 
every other subcarrier in the frequency domain has to be 
set to zero, which is a known Fourier transform property. 

SUMMARY 

[0010] An IEEE 802.16a OFDM system, which uses two 
identical training sequences of one preamble over one OFDM 
symbol duration, may require more sophisticated channel 
estimation techniques at the receiver, compared to a 
system, such as IEEE 802.11a/g OFDM, where training 
sequences of two OFDM symbol durations are transmitted per 
packet. Sophisticated channel estimators, such as the 
Minimum Mean- Square Error (MMSE) and Least -Square (LS) 
estimators have been proposed. However, these estimators 
may be complex to implement because they require many 
multiplication operations. Since training symbols are sent 
for alternate subcarriers {i.e., P k = 0 in equation (l) for 
alternate subcarriers) , interpolation techniques are used 
to estimate the channel response of these zeroed out 
subcarriers . 

[0011] An "upsampling" technique may interpolate the 
channel response. This technique, however, requires 
multiple Fast Fourier transforms (FFT) computations. 
Furthermore, since typical OFDM systems such as IEEE 
802.11a/g and IEEE 802.16a have many guard subcarriers in 
an OFDM symbol that are set to zero, this technique may not 
provide optimum performance. 
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[0012] Other proposed approaches use various filtering 
techniques to compute the interpolated points of an 
estimated channel response. These techniques may typically 
require filtering and multiplication operations. 

[0013] The present application relates to a relatively 
low- complexity channel estimation technique for OFDM 
systems. The technique may use only addition operations. 
The technique may combine (a) linear interpolation with (b) 
relatively low or reduced complexity adaptive filtering to 
achieve a desired level of performance. 



BRIEF DESCRIPTION OF THE DRAWINGS 

[0014] Fig. 1A illustrates a wireless communication 
system with a first transceiver and a second transceiver. 

[0015] Fig. IB illustrates an example of a data 
burst/packet transmitted by the transmitter to the receiver 
in Fig. 1A. 

[0016] Fig. 1C illustrates a set of subcarriers for a 
preamble "p" at the transmitter of Fig. 1A and how the 
subcarriers are transformed and transmitted and how the 
subcarrier channels are estimated. 

[0017] Fig. 2 illustrates a transmit portion of the 
first transceiver and a receive portion of the second 
transceiver in Fig. 1A. 

[0018] Fig. 3 illustrates a technique of using the 
system of Fig. 2. 
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[0019] Fig. 4 shows a table of number of updates and 
adaptation coefficient settings to approximate an optimal 
Kalman filter, which may be used by the system of Fig. 2. 

DETAILED DESCRIPTION 

[0020] Fig. 1A illustrates a wireless communication 
system 130 with a first transceiver 100 and a second 
transceiver 102. The first transceiver 100 may be 
designated a "transmitter" because it first transmits 
signals to the second transceiver 102, which may be 
designated a "receiver." However, both transceivers 100, 
102 may transmit and receive wireless signals, as shown by 
the transmit portions 101A, 101B and receive portions 103A, 
103B. Each transceiver 100, 102 may have one or more 
antennas 104, 106. 

[0021] Fig. IB illustrates an example of a data 
burst/packet 120 transmitted by the transmitter 100 to the 
receiver 102 in Fig. 1A. 

[0022] The transmitter 102 may send a "preamble" to the 
receiver 102. A preamble is an initial set of symbols sent 
by a transmitter for establishing frequency synchronization 
and channel estimation at a receiver. A "symbol" in 
digital transmission is a recognizable electrical state, 
which is associated with a "signal element" {an electrical 
signal within a defined period of time) . For example, a. 
signal element in binary transmission is represented as one 
of two possible states or symbols: 1 or 0 . 
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[0023] Fig. 1C illustrates a set of subcarriers for a 
preamble "p" at the transmitter 100 and how the subcarriers 
are transformed, transmitted and the subcarrier channels 
are estimated at the receiver 102. A subcarrier (also 
referred to as "carrier" or "tone") is a frequency domain 
unit, which may contain a data, pilot, null or training 
symbol (i.e., a symbol with a sequence of known transmit 
signals used for training at the receiver) . 

[0024] Fig. 2 illustrates a transmit portion 101A in the 
transmitter 100 and a receive portion 103B in the receiver 
102 in Fig. 1A. The transmit portion 101A includes a 
multiplexer 206, a modulator 208, an N-IFFT (N-point 
Inverse Fast Fourier Transform) module 210, and a cyclic 
prefix insertion module 212. The multiplexer 206 receives 
a preamble 204 and a "data in" stream, which may be 
scrambled, coded and interleaved. 

[0025] In an IEEE 802.16a OFDM (Orthogonal Frequency 
Division Modulation) system, the preamble P 2 04 contains +1 
or -1 in the 100 subcarriers corresponding to the 100 
subcarriers that would normally carry data and pilot 
symbols during regular transmission. The N-IFFT module 210 
performs a N-point Inverse Fast Fourier Transform. "N" 
represents a number of subcarriers and a number of points 
in the IFFT (or FFT) . The cyclic prefix insertion module 
212 inserts a cyclic prefix (CP) . The transmit portion 
101A transmits a signal 214 with the preamble across a 
wireless channel 155 to the receive portion 103B. 
Modulation by a carrier frequency is not shown for 
simplicity. 
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[0026] The receive portion 103B includes a channel 
estimator 216, an N/2-FFT {Fast Fourier Transform) module 
220, a decision slicer 226, a demodulator 228, an N-FFT 
module 230, and a cyclic prefix removal module 232. The 
cyclic prefix removal module 232 receives a transmitted 
preamble and later receives data from the transmit portion 
101A. The cyclic prefix removal module 232 removes a 
cyclic prefix, and outputs a received preamble sequence P 
231 to the N/2-FFT module 220 in a training mode. In a 
data mode, the cyclic prefix removal module 232 removes a 
cyclic prefix and uses a switch 233 to output data y from a 
received data packet to the N-FFT module 230. 

[0027] The N/2-FFT and N-FFT modules 220, 230 perform 
128-point and 256-point fast Fourier transforms {when N = 
256) , respectively. The N/2-FFT module 220 outputs a 
received preamble P 217 to the channel estimator 216. 
The N-FFT module 23 0 outputs a data symbol Y to the 
demodulator 228. The demodulator 228 outputs a demodulated 
data symbol X 227 to the channel estimator 216 and the 
decision slicer 226. The decision slicer 226 may be any 
type of decision maker, e.g., hard decision or decode 
decision. 

[0028] After N/2 FFT transformation, the channel 
estimator 216 may recognize that +/-1 preamble P 204 was 
sent by the transmit portion 101A and may use the pre- 
stored preamble 2 04 to derive a channel estimate for each 
subcarrier. The channel estimator 216 has access to the 
preamble P 204, which is pre-stored in the receive portion 
103B, as shown in Fig. 2. The channel estimator 216 
compares the pre-stored preamble P 204 to the received and 
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processed P 217. The pre-stored preamble P 204 and the 
received preamble P 217 provide initial channel estimates 
for the corresponding subcarrier. The channel estimator 
216 outputs an estimated channel frequency response H 222 
for the corresponding subcarriers. Interpolation and other 
techniques can be used to obtain channel estimates for the 
subcarriers that do not have a training symbol sent in the 
preamble as described later in the text. 

[0029] The decision slicer 226 may be used to refine 
channel estimates, as described below. The decision slicer 
226 may output decoded data dec{X} 224 to the channel 
estimator 216, as well as a deinter leaver, a channel 
decoder and other components 236. 

[0030] A "training sequence" is a sequence of known 
transmit data in the time domain for establishing 
communications between a transmitter 100 and receiver 102 
for a channel. To perform frequency synchronization in an 
IEEE 802.11a OFDM system, the transmit portion 101A 
transmits two repeated, identical, time domain training 
sequences in two preambles P 204 over two OFDM symbol 
durations. An OFDM symbol (Figs. 1B-1C) in the frequency 
domain is made up of data subcarriers, pilot subcarriers 
and null subcarriers, and an OFDM symbol duration is a 
length of an OFDM symbol in the time domain (Fig. 1C) . The 
number of subcarriers (Fig. 1C) determines an FFT size. The 
receive portion 103B computes the correlation between the 
two received time domain training sequences to obtain 
timing information, correct frequency offset, etc. The 
IEEE 802.16a uses one OFDM symbol duration to transmit the 
preambles in order to obtain a repeated sequence in the 
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time domain with an OFDM symbol period. Training symbols 
are sent for "even" subcarriers only. The OFDM symbol is 
further described on pages 145-147 of the IEEE 802.16a 
draft amendment. 



can be used for channel estimation may have a magnitude of 
l and pseudo random phases for "even" subcarriers. In the 
frequency domain, the preamble 204 may be expressed as: 



[0032] As shown, only one hundred of the two hundred 
non- guard subcarriers (of the total 256 subcarriers per 
OFDM symbol) are actually excited to be +/- 1. Under IEEE 
802.11a and 802.16a standards, the transmit portion 101A 
zeroes out the guard subcarriers, e.g., subcarriers k = - 
128, -127 to -101 and 101, 102 to 127 for an 802.16a 
system, in an OFDM symbol for pulse shaping and other 
purposes . 

[0033] A problem is the channel estimator 216 only has 
information for the excited subcarriers (every other 
subcarrier) but needs to somehow derive the channel 
estimate for the zeroed out subcarriers (k = + /-l, +/- 
3 # ...+/-99) where a training symbol is not transmitted (for 
an IEEE 802.16a system, the channel estimates in the +/- 
101, + /-102, ... +/-127, -128 subcarriers do not need to 
be computed because they are guard subcarriers and no data 
symbols are carried in these guard subcarriers) . 

[0034] An upsampling approach has been proposed to 
estimate the channel response of these P k = 0 subcarriers. 



[0031] 



The portion of the preamble 204 in Fig. 2 that 



(2) 



I" ±l/o>* = ±2,±4,...±100 

10 fork = 0,±1,±3,...±99,±101,±102,...,±127 -128 
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The channel estimator 216 uses the determined preamble P 
217 to obtain channel estimates for the even subcarriers. 
An N-point IFFT is done with these channel estimates back 
to the time domain. The result is periodic with a period 
of N/2. The last N/2 samples are zeroed out, and a N-point 
FFT is performed to obtain all the interpolated points. 
Interpolation with upsampling may not be accurate because 
the zeroed out subcarriers in the time domain, especially- 
zeroed out subcarriers +/- 101 to +/-127, -128, is similar 
to adding high frequency components. This may degrade 
performance. Upsampling is also complex because it 
involves multiple FFT and IFFT computations. 

[0035] In an embodiment, a relatively low- complexity 
channel estimation technique may be used to derive an 
initial channel estimate with even subcarriers and use 
linear interpolation and reduced complexity adaptive 
filtering to refine the initial channel estimate. 

[0036] Since the transmit portion 101A excites almost 
every other subcarrier, the receive portion 103B may assume 
that the two even subcarriers are correlated and perform a 
linear interpolation. The receive portion 103B examines 
every other subcarrier (e.g., even subcarriers up to +/- 
100} and estimates the intermediate (odd) value. 

[0037] Figure 3 is a flowchart describing a channel 
estimate technique according to an embodiment. The N/2-FFT 
module 220 averages the two received, identical, repeated 
training sequences in the preamble p 231 in the time 
domain by adding the time domain samples and dividing by 
two (block 305) . This may be expressed as: 
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(3) P* = 




for k=0, 1, 2, 



.. * - i. 

2 



where p k is one received preamble sequence 231 in the 
time domain, and pt + Jv/2 is the other received preamble 
sequence 231 in the time domain in Fig. 2. N represents 
the number of subcarriers (e.g., 256) that carry the 
preamble. Because of additive white Gaussian noise (AWGN) , 
an average may be derived to increase the signal-to-noise 
ratio (SNR) of the received training sequence. 

[0038] The N/2-FFT module 220 computes a 128-point fast 
Fourier transform (FFT) of pto yield P* 217 for k = 0, + 
/-l, +/-2,..., +/-63, -64 (up to only 64 because the 
previous operation divided 128 by 2) (block 310) . 

[0039] The N-FFT module 230 computes a 256-point fast 
Fourier transform (FFT) of the data portion y of the packet 
to yield Y . The demodulator 228 demodulates the output of 
the N-FFT module 230 to yield X. 

[0040] The channel estimator 216 computes initial 
channel estimates in the frequency domain for even 
subcarriers up to one hundred (block 315) . The initial 
channel estimates may be denoted as Ha 7 o 222, with '0' 
indicating a training period (i.e., a time instance when 
training occurs) at time 0. A '1* indicates a next time 
instance when data starts. A '2' indicates a next time 
instance when another data packet starts. In this case, 
training symbols are transmitted in the even subcarriers of 
the preamble P k 2 04, as described in equation (2) for the 
IEEE 802.16a system example. 
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(4) K2k,o = — for k=±l, ±2, ±3, ±63, -64, and P 2k * 

P2i 

o. 

where 2k indicates that only even subcarriers can be 
estimated since training symbols are only sent over even 
subcarriers. The value for P*,o is from the channel 155 
received by the receive portion 103B, and the value for 2 04 
is pre-stored at the receive portion 103B. Since P 2 k e {+/- 
1} typically (for example Equation (2)), no actual division 
operations may be needed. 

[0041] For the odd subcarriers, the channel estimator 
216 linearly interpolates two derived adjacent channel 
estimates by adding them and dividing by two (block 320) . 

(5 ) H» + w = (H".° + H 2 * + 2,o) for k=0, ±1, ±2 ±63. 

2 

where 2k, 2k+l and 2k+2 indicate subcarrier 2k, 2k+l 
and 2k+2, and H2k+i,o is the interpolated channel estimate. 
H2* + 2,o is derived from a training symbol in a similar way as 

H 2 *,o- 

[0042] For the IEEE 8 02.16a system, H* for k z 101, and 
k < -101 need not be computed because they are guard 
subcarriers, and no data is sent in these subcarriers (see 
Equation (2) ) . 

[0043] Since a data packet typically contains more than 
one OFDM symbol, the channel estimator 216 may make use of 
decoded data symbols dec{X} 224 to refine and update the 
channel estimates H 222 produced by the channel estimator 
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216. The channel estimator 216 may use adaptive filtering 
techniques, such as the Kalman filter. Least Mean Square 
{LMS) , and exponential update, to update the channel 
estimates H 222 based on many observations (block 325) . 



decision slicer 226 of received decoded data in a data 
packet to continuously update the channel estimate and thus 
improve performance. 



Kalman filter to update or refine the channel estimates. 
The channel update may be expressed as: 



where H*, B is the original channel estimate, "n" is a 
time index and indicates the nth time that channel 
estimates are refined/updated, (1 - l/n) is an assigned 
weight, and l/n is another assigned weight. As time 
increases, i.e., n = 1, 2, 3, 4, etc., more weight is 
placed on Ha, b in equation (6) . Equation (6) takes an 
average. Yk,» is the received value or symbol for a 
particular kth received OFDM subcarrier. At each time n, 
the channel estimator 216 performs channel estimates for 
all k subcarriers that carry data (not guard and DC 
subcarriers, as described above) . The decision slicer 226 
receives X*,» and may make a hard decision to decide which 
constellation point is closest to the received symbol X*,« . 
Dec(Xt,«) is decoded data output from the decision slicer 



[0044] 



The channel estimator 216 may use a decision by 



[0045] 



The channel estimator 216 may use an optimal 



(6) 
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226. Dec(X*,«) is the closest QAM constellation point to 
Xt,« from decision decoding by the decision slicer 226. Xm,* 
may be expressed as: 

(7) Xm = Y*,„/H*.„ 

[0046] Thus, the channel estimator 216 makes an updated 
channel estimate based on the hard decision of the received 



[0047] Equation (6) may be used for interpolated 
subcarriers as well as un- interpolated subcarriers . 
Equation (6) may be used for IEEE 802.11a or 802.16a. More 
data symbols may further refine equation (6) . 

[0048] Instead of using a hard decision output, the 
channel estimator 216 may tolerate a delay and use a Dec ( X ) 
output of a Viterbi decoder 236 to obtain channel estimate 
updates based on a sequence estimator instead of symbol-by- 
symbol or decision decoding. 

[0049] The channel estimator 216 may use a least mean 
square (LMS) technique for channel updates, which may be 
computed by: 



where P k * 0 corresponds to the even subcarriers in 
equation (2) , P k = 0 corresponds to the odd subcarriers in 
equation (2) ) , and u and Ui are update coefficients for the 



packet . 




for P k * 0 
/orP k = 0 
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LMS algorithm for the even (trained) and odd (interpolated) 
subcarriers respectively. 



the hard-decision symbols computed from the demodulated 
symbols, i.e., a processed received symbol with the 
modulating carrier frequency and the effects of the channel 
removed. To compute channel updates based on the hard- 
decision symbols, the channel estimator 216 may use: 



where a and a x are adaptation coefficients with values 
between zero and one, and can be determined through 
simulations (e.g., table in Fig. 3). Different adaptation 
coefficients a and ai, and thus different updates, may be 
used for excited subcarriers when P k # 0 and for non-excited 
subcarriers when P k = 0 . The adaptation coefficients a and 
ax may be set to be a power of 2 to avoid using a multiplier 
or divider. This is an approximation to the optimal Kalman 
filter described above with equation (9) . The computation 
in equation (9) is simpler than equation (6) because a 
divider is used for 1/n in equation (6), but not in 
equation (9) . Even so, in a similar approach 1/n in 
equation (6) may be fixed to 2 3 and n = 8. Similarly, a 
may be fixed to 2" 3 in equation (9) . 

[0051] The adaptation coefficients a and cii may also be 
adaptively changed according to the number of updates to 
provide better estimate updates. An optimal changing point 



[0050] 



Another approach for channel updates is based on 



(9) 




16 



Attorney Docket No.: MP0284/13361-048001 

may be stored in a table, as shown in Fig. 4. Fig. 4 shows 
a table of number of updates and adaptation coefficient 
settings to approximate an optimal Kalman filter. For 
example, if the channel estimator 216 is at the 2000th 
update, the update constant to be used is a = 2" 11 . The 
performance difference between the optimal Kalman filter 
and the described method is less than 0.2 dB in the noise 
variance of the channel estimate. An alternative choice is 
to choose a constant a for all the updates depending on the 
length of the packet . 

[0052] The length of the packet is known at the 
beginning of the packet. This information can be used to 
pick an a to be used for all the updates such that some 
metric, such as bit error rate, is minimized. This method 
may perform worse than the previous method, but gives a 
better performance than choosing a fixed a for all packet 
lengths . 

[0053] Since the interpolated channel estimates can 
limit the performance of the receiver, one may choose to 
adapt only the interpolated channel estimates (i.e., 
setting a = 0) to further reduce the complexity. Depending 
on the channel condition and the SNR operating point, the 
optimal adaptation coefficients a and a± can have different 
values. Furthermore, the optimal adaptation coefficients a 
and ax can also be set differently for, for example, 
different modulation schemes, channelizations and packet 
sizes as this information is known at the receiver. The 
proposed channel adaptation update can also be applied to 
channel estimation in OFDM systems, such as the IEEE 
802.11a/g systems, which do not require interpolation. In 
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this case, the adaptation update can be performed after the 
initial channel estimates of all the data subcarriers . 

[0054] In an alternative approach, a sequence estimator 
(symbol -by- symbol) or channel decoding may be used instead 
of hard-decision symbols. 

[0055] Simulating the performance of using the proposed 
channel estimation technique shows that, for example, the 
decoder performance is about 1 dB away from the ideal case 
when the channel is known at the receiver for a rate -1/2 
16-QAM system in SUI {Stanford University Interim) -3 
channels . 

[0056] A number of embodiments have been described. 
Nevertheless, it will be understood that various 
modifications may be made without departing from the spirit 
and scope of the application. For example, the techniques 
above are implemented with IEEE 802.11a and 802.16a systems 
as examples. Other OFDM or non-OFDM systems and standards 
may be used with the techniques described above. The 
techniques may be applied to point -to-multiple -point 
systems and multiple-in-multiple-out (MIMO) systems. 
Accordingly, other embodiments are within the scope of the 
following claims. 
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